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I. INTRODUCTION 



Recent experimental observations [l|, 0, [3[ intensified the interest in B-meson physics as 
the prime area in which the validity of the Standard Model can be put to a test. In order 
to perform such a test, theoretical predictions for fg, Bg and Am s / Arrid are needed. In 
principle, the lattice approach to QCD allows one to calculate these and other quantities of 
interest from the first principles. However, the majority of lattice results to date suffer from 
substantial systematic errors which make comparison with experimental results inconclusive. 
One of the largest systematic errors comes from the procedure which relates lattice operators 
to the corresponding operators in continuum QCD. 

There are two main sources of error in this matching procedure. First, lattice actions 
used in simulations usually have a reduced symmetry group. The lack of symmetry leads to 
the introduction of additional lattice operators which complicates the matching. Second, the 
matching is usually performed at the one-loop level in perturbation theory. More often than 
not, the one-loop contributions turn out to be quite large. Large one-loop contributions 
undermine the perturbative approach to matching and force the introduction of a large 
systematic error in order to account for possible higher-order corrections. 

In this paper we will demonstrate that both of these issues can be addressed by an 
appropriate choice of lattice action. Specifically, we will focus on the combination of the 
static effective field theory [i| for heavy quarks, the domain-wall formulation 0, @] for 
light quarks and an improved gauge action [7] for gluons. This action provides a good 
approximation for heavy mesons such as B, and it has already been employed in numerical 
simulations 0. 

We will obtain the renormalization constants connecting bare lattice operators to con- 
tinuum QCD operators renormalized in MS(NDR) scheme. In particular, we will consider 
heavy-light bilinears and four-fermion operators relevant for B° — B° mixing. The matching 
will be done at the one-loop level in perturbation theory. We will split the matching into 
two parts. First, the continuum QCD operators need to be related to the continuum static 
theory operators. This step has already been carried out, and we will quote the results here. 
Second, the continuum static theory operators need to be matched to their lattice counter- 
parts. This second step has not been performed before for our choice of lattice action, and 
thus is the main subject of this paper. 

This paper is organized as follows. In Sec. [TT] we introduce our choice of lattice action 
and discuss the corresponding Feynman rules. In Sec. II III the symmetries of the action are 
discussed. In Sec. IIVI we quote the results of matching the continuum QCD and continuum 
static theory operators. In Sec. |V] we will match the static theory bilinears to their lattice 
counterparts. In Sec. |Vl]we present the matching of four-fermion operators as well as Bb- 
In Sec. IVHI we discuss the application of mean-field improvement to our results. 

II. LATTICE ACTION AND FEYNMAN RULES 

Our lattice action consists of three parts: 

S = •S'static + 'S'dW + 'S'gauge; (1) 

where Static is the static-quark action describing heavy quarks, Sew is the domain-wall 
fermion action describing light quarks and Sgauge is the pure gauge action. In this and all 
subsequent sections we set the lattice spacing a equal to 1. 
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A. Static-quark action 



The simplest formulation of the static-effective-field-theory action on the lattice is given 
by0 

Static = - U o( X ~ ®M X ~ ( 2 ) 

x 

where Uq(x — 0) is the gauge link in temporal direction between sites x and x — 0. The 
static-quark field h satisfies 70/1 = h. The static-quark propagator in momentum space is 
given by 

The quark-gluon vertices are 

V*(k, k') = -igT a 5 0fl e- t{ko+k '° )/2 , (4) 

V^(k, k') = - l -g 2 {T\ T b }5 ,6 0u e-^ +k '^ 2 , (5) 
with T a being a generator of the SU(iV c ) color group. We note that it is not difficult to 



extend the results of this paper to the case of the modified static-quark action of Ref. [10 
The results of such extension will be reported elsewhere. 



B. Domain-wall fermion action 

For a quark with mass m the domain- wall fermion action 0, is given by 

N 

Sdw = Yl ^2^s(x)D^( x ,y)ij s ,{y) + ^mg(x)g(x), (6) 

s,s'=l x,y x 

D™(x, y) = D\x, y)5 ss , + D 5 (s } s')5 xy + (M 5 - 5)5 ss ,5 xy , (7) 
D\x, y) = ^ i [(1 + 7/.)^.(«) W» + (1 - ^WW-to] , ( 8 ) 

f PflV (s = 1) 

D 5 (s, s') = { Pfl^i,., + P L S s -i, S ' (Ks<N) , (9) 
[ Pl$n-i, s > (s = N) 

where ip s (x) is a (4+l)-dimensional Wilson-style fermion field. The fifth dimension extends 
from 1 to N and is labeled by s. Pl = (1 — 7s)/2 and Pr = (1 + 7 5 )/2 are the projectors for 
the left- and right-handed spinors. The domain-wall height M5 is a parameter of the theory, 
which can be set ^ M5 ^ 2. Finally, q(x) represents the physical four- dimensional quark 
field constructed from the five- dimensional field ip s {x) at s = 1 and N 

q(x) = P R ij 1 (x) + P L ij N (x), (10) 

q(x)=^ 1 (x)P L + ^ N (x)P R . (11) 
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In momentum space the free field domain-wall Dirac operator is given by 

D (p)st = ^2 h ^ sin st + {W + (p) st + m5 sA 5 Njt )P R + {W~(p) st + m5 S}N 5 lit )P L , 



(12) 



where s and t label coordinates in the fifth dimension and 

/-W(p) 1 

W + (p) st = 



\ 



W(p) 1 

-W(p) ) 



(13) 



St 



( -W(p) 

1 -W{p) 

V 



1 -W(p)J 
W(p) = 1 - M 5 + J^i 1 ~ cos ^)- 



(14) 



(15) 



From now on we will focus on m = case and assume that N is large. We will rely on the 
formalism developed in Refs. 12j, ll3j in our perturbative treatment of domain- wall fermions. 
By inverting the Dirac equation one can write the tree-level DWF propagator as 



S(p)st = (-i^smp^u + W su )G* t {p)P R + (-^sinp^ + W+)G^{p)P L , 



where 



.4 



G*(p) = --[(1 - We- a )e- a{2N - s - l) + (1 - We a )e~ a{s+t) } + Ae 



■a\s-t\ 



(16) 



(17) 
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G L st {p) = - We a )e- a{2N - s - t+2) + (1 - We- a )e~ a{s+t - 2) } + Ae'^-^. (18) 

F 

Note that a, A, F and W in the above formulas all depend on the momentum p as follows: 



cosh(a) 



.4 



1 + W 2 {p) +Y,^ 2 p^ 
2\W(p)\ ' 

1 



(19) 



(20) 



2Wsmh(a) 

F = l-e a W(p). (21) 

The form of the propagator in Eq. (I16p is valid for positive W. For < M 5 < 1, W is 
positive for any p. For 1 < M 5 < 2, W is negative in the small-momentum region. For that 
momentum region we need to adjust the formulas using the substitution 



e— -> -e—. (22) 
Using Eq. ffTBT) and Eqs. fllOllllj) one can obtain the formula for the physical quark propagator 

r^sinp^ 



S q (p) = (q(-p)q(p)} 



(23) 



4 



For one-loop calculations we will also need propagators connecting the physical quark field 
and the domain-wall fermion field 

(q(-p)^M) = i 2j^PE {e -<N-s)p R + e -*(s-Dp L) _ e - a{e -*(N-s )pL + e -«(^i)p H ) f (24 ) 

(M-pMp)) = ^ s ™ p » ( e - aiN - s) P R + e~ a{s - l) P L ) - e- a {e~ a{N - s) P R + e - a ^P L ). (25) 
F 

The relevant quark-gluon interaction vertices are the same as in the N-flavor Wilson case 
with r = — 1: 

(p + p% . . . (p + 



V£(p, P% = -^T a f 7 „ cos v " 2 " /M + ^ sin vfl^i 1 ^, (26) 
V#(P,P% = ^{T a ,T*} (^sin&±^ + cos^±^) (27) 
C. Gauge action 

In this paper we consider the following class of gauge actions: 

-W = -\ ( (1 - 8ci) ^ ReTr[Up] + Cl J2 ReTr[[/ fl ] ) , (28) 
* \ p i? / 

where go denotes the bare lattice coupling, Up is the path-ordered product of links around 
the lxl plaquette P and Ur is the path-ordered product of links around the 1x2 rectangle 
R. For ci = 0, Eq. fl28|) reduces to the standard Wilson plaquette action. In the case 
ci = -1.40686, ci_= -0.331 and c x = -1/12 we have the doubly-blocked Wilson (DBW2) 
[HI, 15], Iwasaki [l6| and Symanzik @, [13] actions respectively. 



For the Wilson plaquette action, the gluon propagator in momentum space is straight- 
forward 

= , v % M = 5 f- (29) 



For improved gauge action the gluon propagator is given by [Tj 

nCl , , \ _ (1 ~ A^(fc))4 sin(^) sin(^) + ^ 4 sin 2 (4f) A^) 

where A C j} {k) is a function symmetric in v and p whose explicit dependence on c\ and is 
given in Appendix. We assumed the Feynman gauge for both propagators. 

III. SYMMETRIES OF THE ACTION 

We will now discuss the symmetries of our action and their consequences for heavy-light 



operators. Similar discussions for different choices of action were presented in Refs. |18l . Il9 . 
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We will use the following notation to denote lattice heavy-light bilinears 



O l f = hTq 



( s 


(T 


= i) 


p 


(T 


= 75) 


< y, 


(r 


= 7/*) 




(r 


= 7m7s) 


T 


(r 


= <v) 



(31) 



A. Chiral symmetry 



The five-dimensional domain- wall fermion field ip s (x) can be subjected to the following 
chiral rotations: 

1> a {x) - e iQ "*> s (x), 1 < * < iV/2, (32) 

ijj s {x) -> e ia «*> fl (x), iV/2 + lO^iV, (33) 

where {t a } are the generators of the SU(iV/) flavor group and iV is the extent of the fifth 
dimension. Together with the definition of the physical quark field q in Eq. (ITUj) . the trans- 
formations in Eqs. (13211331) act as standard chiral symmetry transformations on q. For finite 



iV this symmetry is not exact. However, as was shown in Refs. the degree of the 

symmetry breaking decreases rapidly with increasing N. Since we take N — > 00 limit in 
this calculation, we can consider the symmetry breaking to be negligible. Thus, under chiral 
symmetry transformations we have S — ► iP, — > iy4 M and vice versa. This implies that at 
the level of the static theory-to- lattice matching we should expect Zs = Zp and Zy = Za, 
where Z's are defined in Eq. fl85|) . 



B. Heavy quark symmetries 

First, we note that since the static-quark field h satisfies the field equation joh = h we 
have Vq = S and A = P as well as T j = Vj and T^- = e^Ak. Second, the action in Eq. ([2]) 
as well as the continuum static-effective-field-theory action are invariant under the following 
SU(2) transformation 

h -> e^+wnh, (34) 

where <pj is a parameter and o^i = |[7fc, 7«]- This is the heavy-quark spin symmetry [2I. [23l|. 
As its consequence we have the following transformations S iAj, P — > zV,- and vice versa. 
This implies that at the level of the static theory-to-lattice matching we should expect 
Zs = Za and Zp = Zy. 



IV. MATCHING FULL-THEORY AND STATIC-THEORY OPERATORS 



In this section we will briefly review the relation of full QCD operators renormalized in 
MS(NDR) scheme to the corresponding static theory operators. 
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A. Heavy-light bilinears 



The QCD bilinear Or = bTq, renormalized in MS(NDR) at the scale can be related 
to the corresponding static theory bilinear Or = hTq, renormalized at the scale /i, as follows 



24,123 



Or (/•*&) = Co r (lJ,b, fj)Or(fj) + 0(A QC d/ (35) 

where Co T (Hb, /•*) is a perturbative coefficient encoding the physics between the scales fib and 
[i. In lattice calculations \x naturally assumes the value of the inverse lattice spacing and 
fib is chosen to match the b quark mass. The coefficients Co T for each Or were calculated 
perturbatively at one-loop level in Ref. [I| and at two-loop level in Ref. 1261. The coefficients 
Co T are different for different bilinears. 

B. Four-fermion operator 

The QCD operator relevant for B° — B° mixing, here renormalized according to the 
MS(NDR) scheme, is 

{V . A)(V . A) = [h»(l - 7b)<7] [M 1 - 75)?] ■ (36) 
We need to consider only the parity conserving part of it: 

Ow+aa = (h^q) (h»q) + (h^g) (h^q) • (37) 

The relation to the static theory operators is as follows: 

Ovv+AA(fJ>b) = Zx(fib, v)0 V v+aa(v) + Z 2 (fib, ^)Oss+pp(v) + 0(Aqcd/W), (38) 

where 

Ow+aa = 2 (h^-fq) (h^q) + 2 (h^-f ls q) (h^lso) , (39) 
Oss+pp = 2 (h^q) (h^q) + 2 (h i+) l 5 q) {h H Kq) , (40) 

with 

hF>(x) = e ±imv - x ^ib{x). (41) 



The coefficients were calculated at one-loop in perturbation theory in Refs. [271 . l28l . |29| 
We note that Z 2 is an order g 2 coefficient. 



V. MATCHING STATIC-THEORY AND LATTICE HEAVY-LIGHT BILINEARS 

In this section we will match the MS(NDR) renormalized static theory heavy-light bilin- 
ears introduced in Sec. IIVI to their lattice counterparts. We will use the lattice Feynman 
rules from Sec. [TTJ 

The Feynman graphs representing one-loop corrections to the tree-level heavy-light bilin- 
ears are depicted in Fig. [TJ The light-quark wave function renormalization factor coming 



from the diagrams (a) and (b) has been calculated in Refs. [13j, |30j. The vertex correction 
(c) has not been calculated before for the case of DWF-static quarks. We will present this 
calculation below. The heavy-quark wave function renormalization factor coming from the 
diagrams (d) and (e) has been calculated for the Wilson gauge action in Ref. [9|. We will 
extend the result of Ref. ^ to include the case of improved gauge actions. 
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(d) 



(c) 




(b) 



(o) 



FIG. 1: One-loop corrections to heavy-light bilinears. 



A. Vertex correction 



At the tree level, the Green's function ((h(x)Tq(x))h(y)q(z)) for small external momenta 
p is given by 



((hTq)hq) t 



1 r (l-^o 2 ) 
ip if 1 



(42) 



where 



Wo = W(0) = 1 - M 5 



(43) 

We need to calculate the one-loop correction to this Green's function coming from the graph 
depicted in Fig. [21 As can be seen from Eq. ( 1251) . for small external momenta the external 
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This implies that the one-loop vertex-corrected Green's function will take the form 

((hTq)hq) = — (r + 5T) [ -— -!», (45) 
ip if 

where 

6T = Jj£m k )K(k)D? p (k)r fe^^^WWWW^ + ^o" 1 ^)} . ( 46 ) 



and 



(47) 

First, we focus on the factor in curly brackets in the integrand of Eq. (j46p . Since the light 
quark-gluon vertex V^ a has 7 P and Wilson parts, it is convenient to consider them separately. 
We denote the part coming from the 7 P term as S x and the part coming from the Wilson 
term as S w . Then 

N 

S x (k) = J2(Q(- k )M k ))( w o~ Sp R + <~ lp L)> ( 48 ) 

s=l 
N 

S w (k) = ^(q(-k)iT s {k))(w»- s P L + w^Pr). (49) 

s=l 

Using Eqs. (124112 51) and performing the sum over s we find 

«*) - (51) 

Combining this result with the Feynman rules presented in Sec. HIl the one-loop correction 
to the vertex in Fig. [2] can be written as 

ST = -g 2 C F j[E i_ 6 ~!t + e D oW T { 5 x(^)7pCos(^) + S w (k)ism(^)} , (52) 

where the second Casimir Cp = (N% — 1)/2N C . Performing the sum over p and dropping 
the terms odd in k we can write it as follows: 

ST = TSV r = -J^C F T(I X + I w ), (53) 

where 

r _ 1C 2 /• ^ 2 (M+4E J ^sin 2 (^)cos 2 (f) + E J (l-^)sin 2 (%) 
x J k (w -e«)(W(k)-e-«)Al ' {bA) 

and 
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The integral I w is finite and can be evaluated numerically. The integral I x is infrared 
divergent. We can write it as follows: 

Ix = % + %, (56) 

where the first term is the finite contribution coming from the difference between the im- 
proved and the Wilson gauge action 

^ = 16tt 2 / {G( Cl ,k)-G(0,k)}, (57) 
Jk 

with G(ci, k) representing the integrand in Eq. (j5"31) . The second term 

r = 16vr 2 f G(0, k) = 167T 2 / — , (58) 

is just the Wilson gauge action case added back in. Now, the integral I™ is infrared divergent. 
We will regularize it by including a gluon mass A and extract the finite part using the method 
described in Ref. |3l|. We have 



1-A; 2 ^ 
k 2 (k 2 + A 2 ) 



I w x =? w x -i w x -l^ 2 I U2 ) U2 , / 2V (59) 



where 



and 



J- = I6vr 2 ^ ^— _ —yyy^ _ e _ Q) ; (61) 

both are finite and can be calculated numerically. 

Then the one-loop correction to the vertex is given by 

2 2 

SV r = ^C F [- ln(AV) - I* -X™+I%- I w ] = ^C F [- ln(A 2 a 2 ) + d}. (62) 

Note that neither the constant d nor the coefficient of the logarithm depend on the matrix 
T at the vertex. The numerical results for I w , T™ and 2?f are tabulated in Tab. fl] The 
values of I^ 1 for different gauge actions can be found in Tab. [TTJ The value of the constant 
d is given in Tab. IIHI These and all other numerical integrals in this paper were calculated 



using the Monte Carlo integration routine VEGAS |32j, with 1000000 sample points. 



B. Heavy quark corrections 

Now, let's focus on the heavy quark propagator corrections. The tadpole graph in Fig. [lb 
is given by the following expression 

S tad bo) = - l ~g 2 C F e-™ J D^{k) = -~c?C F e-^T^ (63) 
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where T Cl is the tadpole integral whose numerical values are 



0.154933(3) Wilson, ci = 
j 0.094759(3) Iwasaki, Cl = -0.331 . , 

Cl ~ * 0.062426(2) DBW2, ci = -1.40686 ' ^ ' 



0.128291(3) Symanzik,ci = -1/12 
The formula for the rising sun diagram in Fig. [TH is given by 

-i(po+feo) 



S rs (p ) = ~9 2 C F / ^(fc-g) * ■ (65) 



The heavy quark wave function renormalization is obtained via 

^-l- ^^^ U (66) 

The calculation of the tadpole contribution to Z\ is trivial. In order to calculate the contri- 
bution from E rs we will use a method similar to the one developed in Ref. We write 

<9£ rs 

— i -^ \ po =o = S cxp + Sd, (67) 

where S exp and E D come from differentiating exp[— i(p + k )} and D^ {k — p) in Eq. fl65|) 
respectively. We get 

S ex P = g'C F jD^{k) i _ e ^ ko+e (68) 

Note that 

— _ = lcotP)--. (69) 

1 - e~ lfc ° 2z v 2 ; 2 v ; 

Then Eq. (1681) can be simplified to 

S cxp = \g 2 C F T^ - \g 2 C F T Cl , (70) 

where 

^=£(§^<°'*>- < 7i) 

The numerical results for are 

0.252721(6) Wilson, ci = 

„( 3 ) i 0.164361(5) Iwasaki, ci = —0.331 , > 

c i ~ ^ 0.093831(4) DBW2,ci = -1.40686 ' ^ > 

0.219018(6) Symanzik,ci = -1/12 

For H D in Eq. fl67|) we obtain 

n 7 7 if&Q 

S D = -^ 2 C F y fe «(^)4sin(^)cos(^) T -^ rT -, (73) 



11 



where ^ 

D${k) = dD °^ k \ (74) 

with k = 2sin(fc /2). We calculated D^{k) numerically and symbolically using Mathemat- 
ica. Using Eq. fl69|) and dropping the term odd in ko we can reduce Eq. f J73|) to 



Z D = -2g 2 C F J k D^(k)co S 2 (^). (75) 

Again, we calculate the difference between the improved and the Wilson gauge action nu- 
merically: 

J- = 16 7 r 2 ^(^'(A;) + i ! )cos 2 (^). (76) 

The results are 

( 1.310153(55) Iwasaki,ci = -0.331 
i£ = ) 3.510284(162) DBW2, ci = -1.40686 . (77) 
{ 0.297720(15) Symanzik, Cl = -1/12 

The contribution of the Wilson part to be added back in is 

^l^^i^li^i^-^ <78) 



where To is the Wilson tadpole integral calculated in Eq. (1641) . Once again, we subtract 
numerically the infrared divergent part using the method of Ref . j.3l| : 

6 = 16tt 2 f f-L - e A = 4.791861(251). (79) 



Ai ^ 

Finally, putting together all one-loop contributions to Z 2 we have 

Z h 2 = 1 + 2-^-C F [- ln(A 2 a 2 ) - 1 - I£ + 9 + 4tt 2 T® - 7r 2 T ], (80) 

or 

Z h 2 = 1 + ln(A 2 a 2 ) + e], (81) 

where 

' 24.480 Wilson, d = 



14.883 Iwasaki, c x = -0.331 
4.914 DBW2, ci = -1.40686 
21.223 Symanzik, ci = -1/12 



(82) 



with a maximum error of ±1 in the last digit. Our result for the Wilson case agrees with 
the result of Ref. Note that the correction to Z 2 is significantly reduced in the case of 
Iwasaki and DBW2 action. 

Now, let us discuss the radiative correction to the mass 



5M = -S tad (p = 0) - S rs (p = 0). (83) 
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We obtain 



19.954 Wilson, ci = 



; W _ 9 r o 2 T (3) _ ^ r I 12.977 Iwasaki, Cl = -0.331 . . 

6M ~ 16^ Cp8n Tci " W^° F X I 7.409 DBW2, Cl = -1.40686 ' (84) 

[ 17.293 Symanzik, Cl = -1/12 

which agrees with the Wilson case result in Refs. BBS). 



C. Relation to MS scheme 

Now, we can calculate Zr relating the continuum static theory bilinear Or, renormalized 
in MS(NDR) , to the lattice bilinear Of at as follows: 

OrGu) = (1 - wlY^Z^Zr^, a)O l r *\a), (85) 

where (1 — w; 2 ) -1 / 2 .^ 1 ^ 2 is a DWF-specific factor, whose origin and numerical values are 
discussed in Refs. 13, and 



Z r (fi, a) = 1 + ^C F (- In ^ + D + \{-2 In % + E) + I(ln ^ + F)) 

<>2 ; C F (-\n\ 2 a 2 + d + l(-2\n\ 2 a 2 + e) + h\n\ 2 a 2 + f)). (86) 



167r 2 v 2 y '2 



The continuum static theory constants D = 1, E = and F = 1/2 were calculated in 
Ref. [H. The values for d and e can be found in Tab. lIIII and Eq. (182j) . The lattice light quark 
renormalization factor f = 1/2 — z?,, where 22 was calculated in Refs. 13|, Eliminating; 
A from Eq. (1861) we obtain 



Z r (u, a) = 1 + -^C F (- ln//V + - - d - - - I ). (87) 
' 16tt 2 v 2 p 4 2 2 ; 



Note that the result for Zr does not depend on V since both .D and d are T independent. 
This is consistent with the observation in Sec. IHII that Zr should be the same for all bilinears 
because of the symmetries of the action. 



VI. MATCHING STATIC-THEORY AND LATTICE FOUR-FERMION OPERA- 
TORS 

In Sec. HV] we discussed the relationship between the full theory operator Ovv+aa and 
static theory operators Ovv+aa and Oss+pp- In this section we will match these static 
theory operators renormalized in MS(NDR) scheme to lattice operators. We will also discuss 
the renormalization of Bb- 
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A. Operators matching 



First, we note that the coefficient Z2 in Eq. ([3B1 is already of the order g 2 . Thus, for the 
purposes of matching at the leading order in g 2 it will suffice to match Oss+pp and Ogg +PP 
at tree level. 



Oss+pp(/J<) = Z s O 



lat / 
SS+PP\ 



(88) 



with 



1 

W 5 



(89) 
to the 

lattice has been worked out for the Wilson case in Refs. [35|, [36J. In our case, the mixing is 



being the DWF-specific factor present even at tree level. The matching of O 



VV+AA 



simplified, since the chiral symmetry prohibits mixing with operators of different chirality 
With the exception of the diagrams depicted in Fig. [31 the four-fermions 



diagrams can be reduced to the diagrams already discussed in Sec. [V] Thus, we have 



O 



VV+AA 



(m) 



(90) 



with 



and 



w, 



2W7-1 



167P 



(4 In n 



V 



Dt = c a e t v. 

3 3 3 3 3 J 3 



(91) 



(92) 



The constants d, e and / were discussed in Sec.lVl The constant v is the light-light correction 



depicted in Fig. [3^, which was calculated in Refs. [13|, [30j, [39J- In the notation of Ref. [13 



v = —Vs t p- The constant c is the heavy-heavy correction depicted in Fig. [3b, which we 
calculate below. 




(a) 



(b) 



FIG. 3: Four-fermion specific one-loop corrections. 



Omitting the trivial T structure, the contribution from the graph in Fig. [3b is given by 
the following expression: 



5V 



-iko 



g 2 C F / D*(k)j- - h ^. 



(93) 
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Noting that 

e -ifco i 

(1 - e- ifc °)2 = ~4sin 2 (fc /2)' (94) 
we integrate Eq. ( l93l) by parts and arrive at 

6V = -2g 2 C F ^^'(fc)cos 2 (^), (95) 

which is identical to Eq. (1751) in Sec. [V] Using our results from Sec. |V]we can immediately 
write 

SV = 2-^-C F [-\n\ 2 a 2 - 1 - I* + Q - tt 2 T ] = t|^C f [-2 In AV + c], (96) 

where 

' 4.525 Wilson, d = 

1.905 Iwasaki, ci = —0.331 , Q „, 
-2.495 DBW2,ci = -1.40686 ' ^ ' 

3.930 Symanzik, ci = — 1/12 

with a maximum error of ±1 in the last digit. Our result for the Wilson case agrees with 



the result of Ref. 35 



B. Bb renormalization 



In QCD, the B- meson parameter -Bb is defined as 



(*10%gl*) 



Bb - ^J^f^- 08) 



Combining the results of this section with the results of Sec. IIVI and Sec. |V] we can relate 
Bb to lattice as 

_ z x z L lat z 2 z s lat 

~~ 7^2 7 2 -^VV+AA + n2 7 2 SS+PPi l yy J 

where 



|K0|^*|S)|2mB 



Note that the Bb renormalization does not depend on (1 — Wq) 1 Z w l since it cancels in the 
ratios. 



VII. MEAN-FIELD IMPROVEMENT 



Finally, let us briefly discuss the application of mean-field improvement to our results. 



We follow the much more detailed discussion in Refs. [13, The main reason for the 

improvement program is the fact that for some values of M 5 , the DWF-specific factor Z w 
becomes quite large This problem can be circumvented by substituting wq and Z w 

with: 



MF 



w Q + 4{l-u), 



;ioi) 
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z w = Z w \ W0=w mf + — MF 2(1 - it), (102) 

where u = P 1 / 4 , with P being the value of the plaquette. 

The mean field improvement will also affect other perturbative constants. The constant 



/ = I/2-22 in Eqs. (1H71I921 should be replaced by / MF = l/2-^ F , with zf ¥ of Refs. |13|,l3Q 



Because of the change in w , the value of the constant d in Eq. (1621) should also change. 
One can obtain the values of d MF from d by observing that d MF (M 5 ) = d(M 5 ), where 
M 5 = M 5 — 4(1 —«). Similarly, the constant v in Eq. dS2} should be replaced by v MF (M ; 



v(M 5 ). For the heavy quark, mean- field improvement has no effect on the wave function 



renormalization as was pointed out in Refs. [40j, |41]. However, one should be mindful of the 



subtlety involving different normalization conditions for Z\ discussed in Refs. 0, [33[ . The 
normalization condition consistent with the mean-field improvement necessitates the use of 
e — 5M instead of e in Eqs. ( l8Tf92l . The constant c in Eq. fl97|) is unaffected by mean-field 
improvement. Finally, the light-quark field should be multiplied by \/u, which leads to the 
multiplication of Z r by y/u and of Zl,s by u - 



VIII. CONCLUSIONS 

In this paper we have calculated one-loop renormalization constants for operators com- 
bining static heavy and domain-wall light quarks. We obtained results for Wilson, Iwasaki, 
DBW2 and Symanzik gauge actions. We have confirmed that all bilinears regularize by a 
single constant when matched to the static theory. We have shown that the combination 
of domain-wall fermions and an improved gauge action reduces the size of perturbative cor- 
rections significantly compared to the Wilson case. This suggests that our choice of action 
is quite suitable for the determination of fs, P>b and Am s /Amd with increased precision. 
Our results provide the necessary connection between those phenomenological quantities 
and lattice calculations. 
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IX. APPENDIX 

The improved gauge propagator in momentum space is given by [1 



(103) 
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where a is the gauge parameter, k^ = 2 sin(A; At /2), k 2 = 4^] sin 2 (/c M /2) and 



+ c2 ((e^ 4 ) +* a E*JE^+(* 2 ) 2 n^)]> ( io4 ) 

P P T^P:^ Pt^Pi^ 



with 



aw = (£ 2 - Cl p - ci ((^ 2 ) 2 +E^) + r 2 (^ 2 ) 3 + 2 E^ - p E **) 

p r r r 

-4<*e*jii?- (ios) 

P T^P 

The different choices for the parameter c\ correspond to different gauge actions discussed in 
the main text of this paper. In the case C\ — 0, the propagator in Eq. (11031) reduces to the 
standard Wilson plaquette action propagator. 
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TABLE I: Results for vertex integrals defined in Sec. El 


M 5 




TW 

x x 


TW 

x x 


0.05 


-4.138229(319) 


-1.909304(325) 


-0.703771(17) 


0.10 


-4.034432(236) 


-2.030270(283) 


-0.712970(18) 


0.15 


-3.948298(192) 


-2.133953(315) 


-0.722356(18) 


0.20 


-3.872189(163) 


-2.228815(258) 


-0.731938(18) 


0.25 


-3.802874(143) 


-2.316957(249) 


-0.741722(18) 


0.30 


-3.738464(129) 


-2.400620(227) 


-0.751716(18) 


0.35 


-3.677754(117) 


-2.480180(264) 


-0.761929(19) 


0.40 


-3.619923(108) 


-2.558862(216) 


-0.772367(19) 


0.45 


-3.564378(100) 


-2.634128(243) 


-0.783041(19) 


0.50 


-3.510664(93) 


-2.709506(214) 


-0.793959(19) 


0.55 


-3.458427(87) 


-2.782736(229) 


-0.805132(20) 


0.60 


-3.407373(81) 


-2.855975(230) 


-0.816570(20) 


0.65 


-3.357255(75) 


-2.928999(222) 


-0.828285(20) 


0.70 


-3.307856(70) 


-3.001320(248) 


-0.840288(20) 


0.75 


-3.258987(65) 


-3.074116(255) 


-0.852592(20) 


0.80 


-3.210479(60) 


-3.147106(260) 


-0.865211(21) 


0.85 


-3.162165(55) 


-3.220438(268) 


-0.878159(21) 


0.90 


-3.113899(50) 


-3.294598(273) 


-0.891453(21) 


0.95 


-3.065532(44) 


-3.369698(280) 


-0.905107(21) 


1.00 


-3.016959(36) 


-3.445750(258) 


-0.919141(22) 


1.05 


-2.967904(40) 


-3.523295(291) 


-0.933573(22) 


1.10 


-2.918434(44) 


-3.602252(251) 


-0.948424(22) 


1.15 


-2.868167(47) 


-3.682678(253) 


-0.963696(24) 


1.20 


-2.816869(48) 


-3.764914(309) 


-0.979474(23) 


1.25 


-2.764631(52) 


-3.849567(315) 


-0.995722(23) 


1.30 


-2.711104(56) 


-3.936509(322) 


-1.012510(29) 


1.35 


-2.656052(60) 


-4.026280(329) 


-1.029829(30) 


1.40 


-2.599188(65) 


-4.119041(336) 


-1.047733(30) 


1.45 


-2.540188(72) 


-4.215449(345) 


-1.066257(31) 


1.50 


-2.478673(79) 


-4.315845(352) 


-1.085443(31) 


1.55 


-2.414181(87) 


-4.421470(280) 


-1.105334(32) 


1.60 


-2.346139(97) 


-4.531634(286) 


-1.125985(35) 


1.65 


-2.273804(108) 


-4.647352(381) 


-1.147436(33) 


1.70 


-2.196232(120) 


-4.770915(392) 


-1.169762(34) 


1.75 


-2.112133(135) 


-4.903247(399) 


-1.193028(34) 


1.80 


-2.019661(155) 


-5.046261(405) 


-1.217309(35) 


1.85 


-1.916047(182) 


-5.202899(419) 


-1.242693(36) 


1.90 


-1.796524(223) 


-5.378317(435) 


-1.269279(37) 


1.95 


-1.651230(292) 


-5.583553(420) 


-1.297141(31) 
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TABLE II: Results for 1^ integral denned in Eq. 1(57)) ■ 



M 5 


Iwasaki 


DBW2 


Symanzik 


0.05 


0.246331(9) 


0.408168(15) 


0.100822(4) 


0.10 


0.249861(9) 


0.414136(16) 


0.102247(4) 


0.15 


0.253496(9) 


0.420275(16) 


0.103714(4) 


0.20 


0.257237(10) 


0.426590(16) 


0.105225(4) 


0.25 


0.261091(10) 


0.433088(16) 


0.106782(4) 


0.30 


0.265062(10) 


0.439774(17) 


0.108386(4) 


0.35 


0.269153(10) 


0.446657(17) 


0.110041(4) 


0.40 


0.273371(10) 


0.453742(17) 


0.111747(4) 


0.45 


0.277720(10) 


0.461039(18) 


0.113507(4) 


0.50 


0.282207(11) 


0.468556(18) 


0.115323(4) 


0.55 


0.286837(11) 


0.476302(18) 


0.117199(5) 


0.60 


0.291617(11) 


0.484287(18) 


0.119138(5) 


0.65 


0.296554(11) 


0.492522(19) 


0.121141(5) 


0.70 


0.301655(11) 


0.501018(19) 


0.123212(5) 


0.75 


0.306929(12) 


0.509787(19) 


0.125355(5) 


0.80 


0.312383(12) 


0.518844(20) 


0.127573(5) 


0.85 


0.318029(12) 


0.528201(20) 


0.129871(5) 


0.90 


0.323874(12) 


0.537875(21) 


0.132253(5) 


0.95 


0.329931(13) 


0.547883(21) 


0.134722(5) 


1.00 


0.336211(13) 


0.558242(21) 


0.137285(5) 


1.05 


0.342726(13) 


0.568972(22) 


0.139946(6) 


1.10 


0.349490(13) 


0.580095(22) 


0.142712(6) 


1.15 


0.356519(14) 


0.591634(23) 


0.145589(6) 


1.20 


0.363828(14) 


0.603614(23) 


0.148583(6) 


1.25 


0.371435(14) 


0.616064(24) 


0.151703(6) 


1.30 


0.379360(15) 


0.629014(24) 


0.154956(6) 


1.35 


0.387625(15) 


0.642497(25) 


0.158353(6) 


1.40 


0.396252(16) 


0.656551(25) 


0.161902(7) 


1.45 


0.405268(16) 


0.671216(26) 


0.165615(7) 


1.50 


0.414702(16) 


0.686539(27) 


0.169504(7) 


1.55 


0.424586(17) 


0.702570(27) 


0.173584(7) 


1.60 


0.434956(17) 


0.719365(28) 


0.177868(8) 


1.65 


0.445852(18) 


0.736989(29) 


0.182376(8) 


1.70 


0.457321(19) 


0.755513(30) 


0.187125(8) 


1.75 


0.469413(19) 


0.775019(31) 


0.192138(8) 


1.80 


0.482186(20) 


0.795598(32) 


0.197440(9) 


1.85 


0.495709(21) 


0.817357(33) 


0.203059(9) 


1.90 


0.510059(22) 


0.840417(35) 


0.209029(10) 


1.95 


0.525326(24) 


0.864919(37) 


0.215387(10) 
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TABLE III: Results for d defined in Eq. ([62]). 



M 5 


Wilson 




Iwasaki 


DBW2 




Symanzik 


0.05 


5.343762 


[661) 


5 


097431 


;670) 


4.935594 


;676) 


5 


242940 


;665) 


0.10 


5.351732 


;537) 


5 


101871 


^546) 


4.937596 


[553) 


5 


249485 


[541) 


0.15 


5.359895 


;525) 


5 


106399 


[534) 


4.939620 


[541) 


5 


256181 


[529) 


0.20 


5.369066 


;439) 


5 


111829 


[449) 


4.942476 


[455) 


5 


263841 


[443) 


0.25 


5.378109 


;4io) 


5 


117018 


[420) 


4.945021 


[426) 


5 


271327 


[414) 


0.30 


5.387368 


;374) 


5 


122306 


[384) 


4.947594 


[391) 


5 


278982 


[378) 


0.35 


5.396005 


;4oo) 


5 


126852 


[410) 


4.949348 


[417) 


5 


285964 


[404) 


0.40 


5.406418 


;343) 


5 


133047 


[353) 


4.952676 


;360) 


5 


294671 


[347) 


0.45 


5.415465 


;362) 


5 


137745 


;372) 


4.954426 


;380) 


5 


301958 


;366) 


0.50 


5.426211 


;326) 


5 


144004 


;337) 


4.957655 


[344) 


5 


310888 


;330) 


0.55 


5.436031 


;336) 


5 


149194 


;347) 


4.959729 


[354) 


5 


318832 


[341) 


0.60 


5.446778 


[331) 


5 


155161 


;342) 


4.962491 


[349) 


5 


327640 


;336) 


0.65 


5.457969 


[317) 


5 


161415 


;328) 


4.965447 


;336) 


5 


336828 


;322) 


0.70 


5.468888 


;338) 


5 


167233 


[349) 


4.967870 


[357) 


5 


345676 


[343) 


0.75 


5.480511 


[340) 


5 


173582 


[352) 


4.970724 


[359) 


5 


355156 


[345) 


0.80 


5.492374 


[341) 


5 


179991 


[353) 


4.973530 


[361) 


5 


364801 


[346) 


0.85 


5.504444 


[344) 


5 


186415 


[356) 


4.976243 


[364) 


5 


374573 


[349) 


0.90 


5.517044 


[344) 


5 


193170 


[356) 


4.979169 


[365) 


5 


384791 


[349) 


0.95 


5.530123 


[345) 


5 


200192 


[358) 


4.982240 


;366) 


5 


395401 


[350) 


1.00 


5.543568 


;316) 


5 


207357 


;329) 


4.985326 


;337) 


5 


406283 


[321) 


1.05 


5.557626 


[353) 


5 


214900 


;366) 


4.988654 


[375) 


5 


417680 


[359) 


1.10 


5.572262 


[317) 


5 


222772 


^330) 


4.992167 


;339) 


5 


429550 


;323) 


1.15 


5.587149 


[324) 


5 


230630 


^338) 


4.995515 


[347) 


5 


441560 


;330) 


1.20 


5.602309 


;380) 


5 


238481 


[394) 


4.998695 


[403) 


5 


453726 


;386) 


1.25 


5.618476 


;390) 


5 


247041 


;404) 


5.002412 


[414) 


5 


466773 


;396) 


1.30 


5.635103 


;407) 


5 


255743 


[422) 


5.006089 


[431) 


5 


480147 


[413) 


1.35 


5.652503 


[419) 


5 


264878 


[434) 


5.010006 


[444) 


5 


494150 


[425) 


1.40 


5.670496 


[431) 


5 


274244 


[447) 


5.013945 


[456) 


5 


508594 


[438) 


1.45 


5.689380 


;448) 


5 


284112 


[464) 


5.018164 


[474) 


5 


523765 


[455) 


1.50 


5.709075 


[462) 


5 


294373 


;478) 


5.022536 


;489) 


5 


539571 


;469) 


1.55 


5.730317 


;399) 


5 


305731 


[416) 


5.027747 


[426) 


5 


556733 


[406) 


1.60 


5.751788 


[418) 


5 


316832 


[435) 


5.032423 


[446) 


5 


573920 


[426) 


1.65 


5.773720 


[522) 


5 


327868 


[540) 


5.036731 


[551) 


5 


591344 


[530) 


1.70 


5.797385 


[546) 


5 


340064 


[565) 


5.041872 


[576) 


5 


610260 


[554) 


1.75 


5.822352 


[568) 


5 


352939 


[587) 


5.047333 


[599) 


5 


630214 


[576) 


1.80 


5.848613 


[595) 


5 


366427 


[615) 


5.053015 


;627) 


5 


651173 


[604) 


1.85 


5.876253 


;637) 


5 


380544 


[658) 


5.058896 


;670) 


5 


673194 


[646) 


1.90 


5.905562 


[695) 


5 


395503 


[717) 


5.065145 


;730) 


5 


696533 


705) 


1.95 


5.937642 


[743) 


5 


412316 


;767) 


5.072723 


;780) 


5 


722255 


753) 
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